% Plot phase transition

clc;clear;close all;
cc=[0.0549    0.2706    0.6275
    0.9020         0    0.0627
    0.0980    0.5412    0.1843
    0.8941    0.2314    0.5333
    0.9647    0.3294    0.0510
    0.0941    0.5294    0.9098
    0.5686    0.7490    0.1020
    0.9882    0.7098    0.0588
    0.1020    0.5686    0.4745];
load('./data/QPTtranslation.mat');

stairs(Expression1(:,1),Expression1(:,2),'o-','MarkerSize',6,'LineWidth',2,'Color',cc(1,:))
%plot(Expression1(:,1),Expression1(:,2),'MarkerSize',12,'LineWidth',2,'Color',cc(1,:))
xlabel('\delta_x/H_0','LineWidth',4);hold on;
%ylabel('DD invariant / observable','LineWidth',4);
%plot(ones(100),linspace(-0.2,1.2,100),'k--');
%plot(-1*ones(100),linspace(-0.2,1.2,100),'k--')

ylabel('$$\mathcal{DD}$$ invariant','LineWidth',4,'Interpreter','latex');
set(gca,'FontSize',18);
box off;
xlim([-2,2]);ylim([0,1.1]);
